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ABSTRACT 

We report the discovery of two millisecond pulsars in a search for radio pulsations at the 
positions of Fermi Large Area Telescope sources with no previously known counterparts, using 
the Nancay radio telescope. The two millisecond pulsars, PSRs J2017+0603 and J2302+4442, 
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have rotational periods of 2.896 and 5.192 ms and are both in binary systems with low-eccentricity 
orbits and orbital periods of 2.2 and 125.9 days respectively, suggesting long recycling processes. 
Gamma- ray pulsations were subsequently detected for both objects, indicating that they power 
the associated Fermi sources in which they were found. The gamma-ray light curves and spectral 
properties are similar to those of previously-detected gamma-ray millisecond pulsars. Detailed 
modeling of the observed radio and gamma-ray light curves shows that the gamma-ray emission 
seems to originate at high altitudes in their magnetospheres. Additionally, X-ray observations 
revealed the presence of an X-ray source at the position of PSR J2302+4442, consistent with 
thermal emission from a neutron star. These discoveries along with the numerous detections of 
radio-loud millisecond pulsars in gamma rays suggest that many Fermi sources with no known 
counterpart could be unknown millisecond pulsars. 

io Subject headings: pulsars: general — pulsars: individual (J2017+0603, J2302+4442) — gamma 

a rays: general 



Introduction 



Durin g its first year of act ivity, the Large Area Telescope (LAT) aboard the Fermi Gamma- Ray Space 
Telescope (jAtwood et al. 2009f) firmly established millisecond pulsars (MSPs) as bright source s of gamma 



rays, w it h the d etection of pulsed emission from at least nine Galactic disk MSPs above 0.1 GeV (jAbdo et al 



2009dl laL lioiObh . Normal pulsa rs were already establis hed as an important class of gamma-ray sources by 
previous experime nts (see e.g. iThompson et al] Il999l ). The First Fermi Catalog of gamma-ray pulsars 
(jAbdo et al.ll2010ff ) tabulated the properties of 46 pulsars, including eight millisecond pulsars. In addition, 
the LAT has observed gamma-ray emission from seve ral globular clusters (GC s) with spectral properties 
that are consistent with those of populations of MSPs ( Abdo et al. 2009c . 2010af ) and thus of the flux being 
due to the combined MSPs in the cluster. 

Millisecond pulsars are rapidly rotating neutron stars (with rotational period of few tens of millisec- 
onds) with very small spin- down rates (P < 10~ 17 ). They are thought to have acquired their high rota- 
tional rate by accretion of matter, and thereby transfe r of angular momentum, from a binary companion 



Bisn ovatvi-Kogan fc Komberall974l : lAlpar et al.lll982f ). which is now supported by observational evidence 
(lArchibald et al. 2009). About 10% of th e ~ 2000 known pulsars are MSPs, either in the Galactic disk or in 



globular clusters jManchester et ah] l 2005l) . Estimates for the Galactic population of MSPs range from 40000 
to 90000 objects (see iLorimerl 12008. and references therein). A small fraction of these have large enough 
spin-down luminosities E and small enough distances d to be detectable by the LAT. The minimum \f~E~ /d 2 
of pulsars in the Fermi First Pulsar Catalog is 0.1% of the value for Vela. Furthermore, the sparsity of the 
photons recorded by the LAT makes MSPs much easier to discover at rad io wavelengths than in gamma rays 
(for a discussion of blind period searches of gamma-ray pulsars, see e.g. Abdo et al. 2009bl and references 
therein). However, also when blindly searched in the radio band, the MSPs are difficult targets. On one hand 
they are faint sources so that their detection generally requires long exposures with large radio telescopes. 
In addition, most MSPs are in binary systems so the orbital motions need to be taken into account when 
searching for pulsations, introducing additional parameter combinations, and therefore making data analyses 
computationally intensive and searches less sensitive than for normal pulsars. 



Radio emission from pulsars is also affected by pulse scattering induced by the ionized component of 
the interstellar medium, with a characteristic timescale r, oc f~ 4 d 2 where / is the observing frequency and 
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d the pulsar distance (jLorimer fc Kramerl 120051 ) . The short rotational periods of MSPs thus introduce an 
observational bias favoring nearby objects. As a consequence of their proximity and their age, they are more 
widely distributed in Galactic latitude than normal pulsars. 



The Fermi Large Area Telescope First Source Catalog (1FGL) ([Abdo et al.ll2010d) has 1451 sources, 
including 630 which are not clearly associated with counterparts known at other wavelengths. The detec- 
tion of nine radio-loud MSPs in gamma rays strongly suggests that a fraction of high Galactic latitude 
unassociated Fermi sources must be unknown MSPs. Such a source of continuous gamma-ray emission can 
be deeply scanned for pulsations at radio wavelengths, resulting in MSP discoveries, provided their radio 



Kerr et al. 


2011: 


Keith et al. 


2011; 



Roberts et al.ll2011 h 



Most high Galactic latitude gamma-ray sources are blazars and other Active Galactic Nuclei (AGNs). 
Fortunately, distinctive indicators of gamma-ray emission from a pulsar are the shape of the spectral emission 
and t he lack of flux var iability in gamma rays. Gamma-ray pulsars indeed exhibit sharp cutoffs at a few 
GeV (jAbdo et al.ll2010ll ). while blazars are known to emit above 10 GeV with no sharp energ y cutoff (Flat 
Spectrum Radio Quasars are well-described by broken power-law spectra) (jAbdo et al. 2010d ). Also, know n 
gamma-ray pulsars are steady sources, whereas blazars show variations of flux over time (jAbdo et alj|2010d) . 
In this expl o ratory study we limited our source discrimination criterion to spectral shapes. As suggested by 
Story et al.1 ( 20071 ). follow-up radio searches of Fermi sources having hard spectra with cutoffs should yield 
discoveries of new MSPs. Gamma-ray variability will be exploited in future studies. 

In this article, we present the observations of pulsar candidates made at the Nangay radio telescope 
that led to the discovery of the MSPs J2017+0603 and J2302+4442 (Section H}. Following the detections, 
we made radio timing observations at the Nangay, Jodrell Bank and Green Bank telescopes (see Sections 13.11 
and !4.ip . The initial ephemerides for these 2.896 and 5.192 ms pulsars in low-eccentricity orbits around light 
companions allowed us to detect gamma-ray pulsations in the data recorded by the LAT. In Sections 13. 3[ 
4.31 and 15.21 we discuss the gamma-ray properties of the two MSPs, compare the measured light curves and 
spectral properties with those of previously observed gamma-ray MSPs. We finally present results of radio 
and gamma-ray light curve modeling in the context of theoretical models of emission in the magnetosphere 
in Section [Ol 



2. Search observations 

The list of 1FGL catalog sources searched for pulsations with the Nangay radio telescope was constructed 
using the following criteria. The radio search was based on a preliminary list of Fermi LAT sources used 
internally by the instrument team. The selection described here is the same as that used, but was applied to 
the 1FGL catalog and yielded the same targets. We first removed gamma-ray sources associated with known 
objects. Sources below —39° in Declination were rejected, as they are not observable with the telescope. 
Sources with Galactic latitudes |fe| < 3° were excluded, being more likely affected by radio pulse scattering 
and also being less accurat ely localized in gam ma rays because of the intense diffuse gamma-ray background 



at low Galactic latitudes (jAbdo et al.ll2010d ). The Nangay beam has a width at half maximum of 4' in 
Right Ascension; therefore we applied a conservative cut by requiring the semi-major axis of the gamma-ray 
source 95% confidence ellipse to be less than 3'. Finally, we selected objects with spectra deviating from 
simple power laws, i.e., showing evidence for a cutoff, and therefore likely pointing to gamma-ray pulsars. 
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For that we excluded so urces with curvatu re indices below 11.34, the limit at which spectra start departing 
from simple power laws (jAbdo e t al. 2010c ). Details on the determination of positions and curvature indices 
of 1FGL sources can be found in lAbdo et al.l (|2010dY 
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From these selection criteria we obtained a list of six sources. Four of them, 1FGL J0614. 1—3328, 
J1231. 1-1410, J1311. 7-3429 and 1FGL J1942. 7+1033, have been searched for pulsations with the Green 
Bank and Effelsberg radio telescopes, and radio pul sars have been detected in the first t wo sources. The 
results of these searches will be reported elsewhere (| Ransom et al.l 1201 It iBarr et al.l 120111 ) . We carried out 
radio observations at the Nangay radio telescope of the other two sources in this list, 1F GL J2017. 3+0603 and 
J2302. 8+4443, using the m odified Berkeley-Orleans-Nangay (BON) instrumentation ( Theureau et al. 20051; 
Cognard fc Theureaul 120061 ) at 1.4 GHz. Instead of doing the usual coherent dedispersion of the signal, the 
code was modified to get a 512 x 0.25 MHz incoherent filterbank sampled every 32 /is. The very first data 
samples were used to determine an amplitude scaling factor and total intensity is recorded as a 4-bit value. 
Observations were usually one hour long, mainly limited by the fact that Nangay is a meridian telescope. 

Data were searched for a periodic dispersed signal using the PRESTO package |Ransom et al.ll2002h . 
After the standard RFI-excision procedure, a total of 1959 dispersion measure (DM) values up to 1244 pc 
cm were chosen to dedisperse the data. Searches for periodicity were done using the harmonic summing 
method (up to eight harmonics). We also searched the data for single pulses, and did not find any. 

An observation of 1FGL J2302. 8+4443 performed on 2009 November 4 revealed a candidate with a 
period of 5.192 ms and a DM of 13.4 pc cm~ 3 . Confirmation observations scheduled at Nangay and Green 
Bank (at 350 MHz) later firmly established this new millisecond pulsar. A week after that first discovery, 
a second candidate in 1FGL J2017. 3+0603 with a period of 2.896 ms and DM of 23.9 pc cm' 3 was also 
confirmed with subsequent Nangay and Green Bank Telescope observations as well as with old observations 
made at the Arecibo telescope. In both cases, substantial variations of the pulsar rotational period were 
observed, indicating orbital motions, as discussed in Sections 13. II and 14.11 



Integrated radio profiles at 1.4 GHz are presented in Figures[T]and[21 The pulse profile of PSR J2017+0603 
is complex and exhibits at least five components. A sharp peak is observed, making PSR J2017+0603 a 
promising addition to pulsar timing array programs. The radio profile of PSR J2302+4442 is broad, with 
at least four pulsed components, three of which form a first structure whose mid-point is separated by ~ 
0.6 rotation from the fourth component. The mean flux density averaged over all observations for the two 
puls ars was determined u sing a calibrated pulse noise diode fired for 10 seconds before each observation 
fsee iTheureau et al.ll201l[ for a description of radio flux measurements with the Nangay radio telescope). 
PSR J2017+0603 presents a mean flux density at 1.4 GHz of 0.5 ± 0.2 mJy, while PSR J2302+4442 is 
brighter at 1.2 ± 0.4 mJy, both being typical values for millisecond pulsars. 



3. PSR J2017+0603 

3.1. Timing observations 

After the initial discovery of PSR J2017+0603, timing observations we re undertaken at t he Nangay 
radio telescope and the Lovell telescope at the Jodrell Bank Observatory (jHobbs et al.ll2004r ). Nangay 
timing observations were done using two different configurations of the BON instrumentation described 
above. Between MJDs 55142 and 552 28 we used the 512 x 0.25 MHz incoherent filterbank at 1334 MHz, 
and the standard coherent dedispersor (|Cog nard et aklkoogh between MJDs 55232 and 55342. The coherent 
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Fig. 1. — Radio and gamma-ray light curves of PSR J2017+0603. The bottom panel shows an integrated 
radio profile at 1.4 GHz with 2048 bins per rotation, recorded with the Nancay radio telescope, based on 
16.2 hours of coherently dedispersed observations. The top three panels show light curves in different energy 
bands (labeled) for gamma-ray events within 0.8° of the pulsar position, with 50 bins per rotation. Two full 
rotations are shown for clarity. See Section [3.31 for details on the determination of background levels, shown 
by horizontal dashed lines. 
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Fig. 2. — Same as Figure [TJ for PSR J2302+4442. The radio profile is based on 20.9 hours of observation. 
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dedispersion is performed in 4 MHz channels over a total bandwidth of 128 MHz centered at 1408 MHz. 
Eighteen times of arrival (TOAs) were recorded with the interbank BON with a mean uncertainty of the 
TOA determination of 6.3 /is, and 19 TOAs were measured with the coherent dedispersor BON with a 
mean uncertainty of 2.5 fis. In addition, 24 radio TOAs were recorded with the Lovell telescope at 1520 
MHz between MJDs 55218 and 55305, with a mean uncertainty of 17.8 fis. These data were used to derive 
an initial timing solution co vering the first seven months post-discovery, using the TEMP02 pulsar timing 
packag^l (jHobbs et al. 2006 ). The dispersion measure was estimated independently: the data recorded with 
the BON backend of the Nangay telescope were cut in four frequency bands of 32 MHz, centered at 1358, 
1390, 1422 and 1454 MHz. We fitted the multi-frequency dataset with the initial timing solution, where the 
DM was left free. We measured DM = 23.918 ± 0.003 pc cm" 3 . 

For this DM and line-of-si ght, the NE2001 mode l of the Galactic distribution of free electron^ assigns 
a distance of 1.56 ± 0.16 kpc ( Cordes fc Laziol 2002 ). Archival optical and infrared images (POSS-II) and 
radio images (NVSS) show no obvious clouds which might indicate ele ctron over densiti es. The line-of-sight 
intersects the Galaxy's Sagittarius arm at about 2 kpc from the Earth |Reid et al.ll2009l) and at the nominal 
DM distance the pulsar environment is not especially crowded. Nevertheless, density variances not modeled 
in NE2001 could change the distance significantly. 

We phase-folded the data recorded by the Fermi LAT using the initial timing solution, and detected 
pulsed gamma-ray emission with high significance. The gamma-ray light curve and spectral properties of 
the MSP are discussed below. However, we observed gradual phase coherence loss for gamma-ray photon 
dates which were earlier than the ephemeris validity interval, defined by the radio observation time span, 
indicating erroneous parameters in the initial timing solution. To enhance the timing solution and make 
it accurate for the entire time rang e of the LAT dat a used here, we extracted TOAs from the gamma-ray 
data using the method described in iRav et al.l (|201ll) . The LAT data were divided in time intervals where 
the gamma-ray pulsation had a significance of at least 3er. For each time interval, we then measured a 
TOA by cross-correlating the observed gamma-ray light curve and a standard template, derived from the 
fraction of the LAT data covered by the timing solution. The pulsar ephemeris was then optimized with the 
gamma-ray and radio TOAs. This procedure was repeated until phase-coherence was ensured over the whole 
LAT dataset. We eventually extracted a total of 10 gamma-ray TOAs between MJDs 54682 and 55294, with 
a mean uncertainty of 49.1 /is. 

The final timing solution was built using radio and gamma-ray TOAs, fitting for the pulsar position, 
rotational period and first derivative, binary parameters and phase jumps between observatories. The dis- 
persion rneasurej^ue_was held fixed at this stage. The low-eccentricity orbit was described using the ELL1 
model (jLange et al.l 120011 ). We corrected for any underestimation of TOA uncertainties and badness of fit 
by using "error factors" (param eters EFAC in TEMP02) on each set of TOAs, following the method de- 
scribed in lVerbiest et al.l (|2009f ). in order to get a reduced x 2 value as close as possible to unity for the entire 
dataset. We obtained a reduced % 2 value of 1.14. The corresponding timing solution is given in Table [TJ 
The spin-down luminosity and magnetic field at the light cylinder de rived from the measured period and 
period derivative are typical of other gamma-ray MSPs detected so far (|Abdo et al.ll2009aLl2010bh . However, 
with a small period derivative of ~ 8.3 x 10~ 21 and at a distance of 1.56 kpc accordi ng to the NE200 1 
model, PSR J2017+0603 is subject to significant contribution from the Shklovskii effect ( Shklovskii 197C)h . 



making the apparent period derivative greater than the intrinsic one, by 2.43 x 10 21 s 1 Pdfi^, where P is 



1 http: / / sourceforge. net / proj ects / tcmpo2 / 



2 Available at | http : / / rsd-www. nrl . navy . mil/7213/lazio /ne_model / 
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the pulsar rotational period in s, d is the distance in kpc, and \it is the proper motion, in mas yr _1 . This 
effect would reduce the true P and thus reduce the calculated spin-down luminosity and magnetic field at 
the light cylinder. In this study we could not measure any significant proper motion, though it may become 
possible with accumulated radio observations. 

Using the measured binary parameters, projected semi-major axis of the orbit, x, and orbital pe- 
riod, P\,, we calculated the mass function in Table [TJ given by /(m p ,m c ) = (m c sini) 3 /(m p + m c ) 2 — 
(4ir 2 c 3 x 3 )/ (GMqP 2 ) where m p is the pulsar mass, m c is the companion mass and i is the inclination of the 
orbit. Assuming an edge- on orbit (i = 90°) and a p ulsar mass of 1.4 M©, we calculate a lower limit on m c 



of 0.18 M Q . As noted in lLorimer fc Kramer! (|2005| ). the probability of observing a binary system with an 
inclination of less than io for a random distribution of orbital inclinations is 1 — cos(io), therefore a 90% 
confidence upper limit on the companion mass can be derived by assuming an inclination angle i of 26°. 
Doing so gives an upper limit of 0.45 M© for the companion mass of PSR J2017+0603. These mass function 
and range of likely companion mass values indicate that the companion star probably is a He-type white 
dwarf. 



3.2. Optical, UV and X-ray analysis 
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We searched for X-ray and optic al/UV counterparts i n Swift (jGehrels et al.ll2004j ) observations obtained 



from Feb.-Mar. 2009. In an XRT dBurrows et al 
measured an upper limit to the 0.5 



20051) image with 16.4 ks of cumulative exposure, we 



keV count rate of < 1.5 counts ks 1 at the pos ition of PSR J2017+ 0603. 

1 (0.3 - 10 keV) from fevans et al.l (|2007h . and 
an appropriate conversion to our choice of energy range, re sults in a flux limit between 0.5 and 8 keV of 



Adopting a flux conversion of 5 x 10 11 erg cm 2 counts 



< 6 x 10 erg cm s . The UVOT (|Roming et al.ll2005j ) images show a relativ ely bright field sou rce (B 



= 19.8 mag, R.A. = 20:17:22.51, Decl. = +06:03:07.7 with < 0.1" uncertainty, from lMonet et al.ll2003l ). that 
is 3.6" away from the pulsar position, which contaminates the photometry. Moving the aperture sufficiently 
to avoid this source, we estimate optical/UV upper limits for the pulsar to be 80 (V), 47 (B), 17 (U), 7 
(Wl), 5 (M2), and 3 (W2) /iJy. All flux upper limits are at the 3cr confidence level. 



3.3. Gamma-ray analysis 

The gamma-ray data recorded by the LAT were analyzed using the Fermi science tools (STs) v9rl6pf|. 
Using gtselect we selected events recorded between 2008 August 4 and 2010 May 26, with energies above 0.1 
GeV, zenith angles < 105°, and within 20° of the pulsar's position. We furthermore selected events belonging 
to the "Diffuse" class of events under the P6-V3 instrument r esponse function (IRFs), those events having 
the highest probability of being photons ( Atwood et al. 20091) . We finally rejected times when the rocking 



angle of the satellite exceeded 52°, required that the DATA.QUAL and LAT .CONFIG are equal to 1 and 
that the Earth's limb did not infringe upon the Region of Interest (ROI) using gtmktime. Finally, we phase- 
folded gamma-ray events using the pulsar ephemeris given in Tabled] and the Fermi plug- in now distributed 
with the TEMP02 pulsar timing package. 

Figure [1] shows radio and gamma-ray light curves of PSR J2017+0603, for gamma-ray events within 
0.8° of the pulsar. Under this cut, most high-energy photons (energies above 1 GeV) coming from the 



3 http: / / fermi.gsfc.nasa.gov/ssc / data/ analysis / scitools / overview.html 
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pulsar are kept, while the contribution o f background emission, mostly present at lower e nergies, is reduced. 
The bin-independent H-test parameter (|de Jager et al.lll989t Ide Jager fc Busching|l2010r ) has a value of 235, 
corresponding to a pulsation significance well above 10c. As can be seen in Figure [TJ the gamma-ray 
pulse profile comprises two close peaks, offset from the radio emission. The absolute phasing in these light 
curves is such that the maximum of the first Fourier harmonic of the radio profile transferred back into the 
time domain is at phase 0. Under that convention, the maximum of the radio profile is at $ r = 0.123 in 
phase. We fitted the gamma-ray light curve above 0.1 GeV using a two-sided Lorentzian function for the 
asymetrical first peak and a simple Lorentzian function for the second peak above constant background. 
For each peak, the peak position $^ and the Full Width at Half-Maximum FWHM^ are listed in Table 
[2j The Table also lists the values of the radio-to- gamma-ray lag S — $i — $ r , and the gamma-ray peak 
separation A = $ 2 — Quoted uncertainties are statistical. For the radio-to-gamma-ray lag S we quote 
a second error bar, reflecting the uncertainty on the conversion of a TOA recorded at 1.4 GHz to infinite 
frequency, due to the uncertainty on the dispersion measure (DM) value given in Table [TJ With 5 ~ 0.22 
and A ~ 0.29, PSR J20 1 7+0603 follows the correlation between 5 and A expected in outer magnetospheric 
models as pointed out by Romani fc Yadigaroglu ( 19951) and effectively observed for currently known gamma- 
ray pulsars (see Figure 4 of lAbdo et al.ll2010a ). However it is interesting to note that this MSP occupies a 
region of the S - A plot where few gamma-ray pulsars were known. 

The spectral analysis was done by fi tting the region around PSR J2017+0603 using a binned likelihood 



method |Cashlll979t iMattox et"al 



1FGL catalog sources (jAbdo et al 



1996), implemented in the pyLikelihood module of the Fermi STs. All 
2010a ) within 15° from the pulsar as well as additional point sources 



found in an internal LAT source list using 18 months of data were included in the model. Sources were 
modeled with power-law spectra, except for PSR J2017+0603 which was modeled with an exponentially cut 
off power-law, of the form: 



dN 
~dE 



N 



-r 
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In Equation (JJJ), Nq is a normalization factor, T denotes the photon index, and E c is the cutoff energy of 
the pulsar spectrum. The parameter (3 determines the steepness of the exponential cutoff. Fermi LAT pulsar 
spectra are generally well-described by a simple exponential model, f3 = 1. The Galactic diffuse emission was 
modeled using the glLiem_v02 mapcube, while the extragalactic diffuse and residual instrument background 
components were modeled using the isotropic-iem_v02 templat^j. Normalization factors and indices for all 
point sources within 7° from PSR J2017+0603 and normalization factors for diffuse components were left 
free. The best-fit values for the photon index and cutoff energy of PSR J2017+0603 for a simple exponentially 
cut off power-law (/3 = 1) are listed in Tabled and the corresponding gamma-ray energy spectrum is shown 
in Figure |31 The first errors are statistical, and the second are systematic. These last uncertainties were 
calculated by following the same procedure as above, but using bracketing IRFs for which the effective 
area has been perturbed by ± 10% at 0.1 GeV, ± 5% near 0.5 GeV and ± 20% at 10 GeV with linear 
interpolations in log space between. We also modeled the millisecond pulsar with a power-law fit, (3 = 0, 
and found that the exponentially cut off power-law model {(3 — 1) is preferred at the 9a level. A fit of the 
pulsar's spectrum with the (3 parameter in Equation (JTJ) left free led to (3 = 1.5 ±0.6. This value is consistent 
with 1 within statistical errors, and the extra free parameter did not improve the quality of the fit, as can be 



'The diffuse models are available through the Fermi Science Support Center (FSSC) (see http://fermi.gsfc.nasa.gov/ssc/) 
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seen in Figure [3j We therefore conclude that the simple exponentially cut off power-law model (with /5 = 1) 
reproduces the present data well. 

With the full spectral model obtained with this analysis and the Fermi ST gtsrcprob, we calculated 
probabilities that each photon originates from the different gamma-ray sources in the ROI. If we denote uji 
as the probability that a given photon has been emitted by PSR J2017+0603, and therefore (1 — U)i) the 
probability that the photon is due to background, then the background level in the considered ROI can 
be estimated by calculating b — J^f^l — where N is the number of photons in the ROI. The back- 
ground levels shown in Figure [T] were calculated with this method, which is more powerful at discriminating 
background events than methods involving surrounding annuli. 

The photon index T and cutoff energy E r measured in this analysis are reminiscent of those of previously- 
detected gamma-ray MSPs (jAbdo et al.ll2009al l2010bh . Integrating Equation ([I]) above 0.1 GeV yields the 
photon flux F and energy flux G given in Table [21 The 1FGL Catalog quotes an energy flux above 0.1 
GeV for 1FGL J2017.3+0603 of (4.5 ± 0.5) xlO" 11 erg cm" 2 s"\ consisten t with the value measured for 
PSR J 2017+0603. Nevertheless, the high-redshift blazar CLASS J2017+0603 (|Mvers et al.ll2003l: lAbdo et al 



2010eh located 2.3' from the pulsar could also contribute to the gamma-ray flux of the 1FGL source. We 
checked that hypothesis by selecting the off-peak region of the spectrum (pulse phases between 0.25 and 
0.75) and by performing a likelihood analysis of the selected data, where the blazar was modeled by a power- 
law. Following this procedure we did not detect any significant emission from the blazar. PSR J2017+0603 
therefore is the natural counterpart of 1FGL J2017. 3+0603. 



256 4. PSR J2302+4442 

257 4.1. Timing observations 

258 Radio timing observations of the pulsar in 1FGL J2302. 8+4443 were conducted at the Nangay radio 

259 telescope in the two configurations described in Section 13.11 the Green Bank Telescope in West Virginia 

260 with the GUPPI backend_|, and the Lovell telescope at the Jodrell Bank Observatory. Between MJDs 55139 

261 and 55218, 29 TOAs were recorded with the fflterbank BON with a mean uncertainty on the determination 

262 of arrival times of 7.6 fxs, while the coherent dedispersor was used to measure 22 TOAs between MJDs 

263 55150 and 55342, with a mean uncertainty of 2.1 /is. The Green Bank Telescope recorded 32 TOAs in two 

264 observation sessions, at MJDs 55095 and 55157, with a mean uncertainty of 5.1 /zs. The Lovell telescope 

265 recorded a total of 38 TOAs at 1520 MHz between MJDs 55217 and 55304, with a mean uncertainty of 

266 20.4 /is. An initial timing solution was built using these radio timing observations and the TEMP02 pulsar 

267 timing package. As with PSR J2017+0603, data recorded with the BON backend were cut in four frequency 

268 bands of 32 MHz, and the multi- frequency TOAs extracted from these observations were used to determine 

269 the dispersion measure. 

270 We measured DM = 13.762 ± 0.006 pc cnT 3 . The NE2001 model assigns this DM and line-of-sight a 
distance of 1.18±g;^ kpc. Again, optical, infrared, and radio images show no clouds. These line-of-sight 

272 and distance place the pulsar within the Orion spur of the Sagittarius arm. As above, unmodeled electron 

273 density variations could change the distance significantly. 

274 We used the initial timing solution to phase-fold the LAT data and detected highly-significant gamma- 



5 https://safe.nrao.edu/wiki/bin/view/CICADA/NGNPP 
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Fig. 3. — Phase-averaged gamma-ray energy spectrum for PSR J2017+0603. The solid black line shows the 
best-fit model from fitting the full energy range with a simple exponentially cutoff power-law functional form 
(/3 = 1). Dashed lines indicate la errors on the latter model. The dot-dashed line represents the spectral fit 
with the (3 parameter left free. Data points are derived from likelihood fits of individual energy bands where 
the pulsar is modeled with a simple power-law form. A 95% confidence level upper limit was calculated for 
any energy band in which the pulsar was not detected above the background with a significance of at least 
2a. 
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ray pulsations. The gamma-ray light curve and spectral properties of the MSP are discussed below. Similarly 
to PSR J2017+0603, we could not fold all LAT data properly using the initial timing solution, as we observed 
loss of phase-coherence for photons recorded before the first radio timing data were taken. Following the 
iterative procedure described in Section 13.11 we extracted TOAs for the gamma-ray data, optimized the 
timing solution by adding the gamma-ray TOAs to the radio dataset, and phase-folded the LAT data until 
we obtained phase-coherence over the entire Fermi dataset described previously. We finally measured nine 
TOAs between MJDs 54682 and 55294 with an uncertainty of 44.6 fj,s. 

The final timing solution obtained by fitting for the pulsar position, rotational period and first time 
derivative and binary parameters is listed in TableQ] The low-eccentricity orbit of PSR J2302+4442 was also 
described using the ELL1 model. The same procedure to correct underestimated TOA uncertainties with 
EFAC parameters as described in l3.1l was used, resulting in a reduced x 2 value of 1.04. Like PSR J2017+0603, 
J2302+4442 is subject to significant contribution from the Shklovskii effect, with a relatively small period 
derivative of ~ 1.33 x 10 -20 . We were not able to measure any significant proper motion with the present 
dataset, however accumulated radio observations may help constrain the Shklovskii contribution. 

Under the assumption of an edge-on orbit and a pulsar mass of 1.4 M©, the lower limit on the companion 
mass is found to be 0.30 M Q . However, assuming an inclination of i = 26° leads to an upper limit of 0.81 M 
for the companion mass, suggesting that the companion star could either be a -ffe-type or a CO-type white 
dwarf. Nevertheless, the orbital period and ecce ntricity of PSR J2302+4442 are in good agreement with the 



Pb - e relationship predicted bv iPhinnevi (|1992h . whereas "intermediate- mass binary pulsars" (IMBPs) with 
heavier companion stars do not necessarily follow the relationship. This suggests that PSR J2302+4442 is in 
orbit with a low-mass -ffe-type companion, and thus that its inclination angle i must be large. Future radio 
timing observations may help determine the companion mass and orbital inclination, via the measuremen t 



of the Shapiro delay (see e.g. lLorimer fc Kramerl 120051 ). As discussed in detail in iFreire fc Wexl (|2010f ) . 



the amplitude of the measurable part of the Shapiro delay for an orbit with medium to high inclination 
is proportional to h 3 = T Q m c x (sin(i)/(l + |cos(i)|)), where T Q = GM Q /c 3 ~ 4.925 490 947 jus. With a 
current average uncertainty on TOAs recorded with the Nangay BON backend of ~ 2.1 fis, we expect the 
Shapiro delay to be measurable for large m c and i values. 



302 4.2. Optical, UV and X-ray analysis 

303 In the Swift XRT image of the PSR J2302+4442 field (9.1 ks summed exposure), there is a marginal 

304 detection (2.6a) of an X-ray source (R.A. = 23:02:47.00, Deck = +44:42:20.7; 90% confidence radius of 6.3") 

305 that is consistent with the pulsar position. The 0.5 - 8 keV flux corresponding to the observed count rate of 

306 (1.0 ± 0.4) counts ks -1 is ~ 4xl0 -14 erg cm~ 2 s _1 (See Section [3~2l for details on the flux conversion). The 

307 optical and UV upper limits at the pulsar position are: 53 (V), 26 (B), 13 (U), 6 (Wl), 4 (M2), and 3 (W2) 
3oa /ijy 

309 On 2009 December 25, while this Fermi LAT source was as yet unidentified, the XMM- Newton satellite 

310 observed the LAT-source field with the EPIC-MOS and -PN cameras in an effort to explore the source region. 
3n We reduced these data with the Science Analysis Software (SAS) version 10.0.0 released on 2010 April 28. 

312 After filtering the observation for intervals of high particle background we were left with good time intervals 

313 consisting of 24.9 ks, 25.1 ks, and 20.8 ks exposures in the EPIC-MOS1, -MOS2, and -PN instruments, 

314 respectively. A number of sources were detected in the field of the Fermi LAT source, as can be seen in 

315 Figure [4] Once the radio pulsar position was refined to the arcsecond level one X-ray source in particular 



- 13 - 



was positionally identified as the likely pulsar candidate and we name this source XMMUJ230247+444219. 

We extracted events from a 50" region around the source from both the MOS1 and MOS2 event files, 
and background events from a 100" region nearby and apparently free of faint X-ray sources but still on the 
same respective MOS CCD chips. For the PN event files, in order to avoid a gap between adjacent CCDs, 
we extracted events from a region only 10" in radius and a background region of radius 80". From the MOS 
instruments we obtain 269 and 262 events, and from the PN we obtain 176 events, respectively, from the 
source regions. This yields, along with the background estimates, a combined detection significance of 13.2cr 
from all three detectors for XMMUJ230247+444219. 



We grouped these events into spectral bins of at least 20 counts per bin and performed a simultaneous 
XSPECH fit to an absorbed power-law model to all three spectra in the 0.4 to 3.0 keV range. This yields 
a power-law index of 5.9 which we regard as unphysical and so we discard thi s model. On the o ther hand, 
an absorbed neutron star hydrogen atmosphere model (phabs x nsatmos, see Heinke et al. 20061) yields an 



acceptable fit, provided that the neutron star mass and radius are fixed at 1.4 Mq and 10.0 km, respectively, 
and the source distance is fixed at the DM value of 1.18 kpc. However, while we obtain an acceptable 
reduced x 2 of 1.032 for 17 degrees of freedom, we measure a column density of Nh = 0.018+q'q* 8 x 10 22 cm" 2 
(90% confidence) meaning that Njj is poorly constrained and consistent with values anywhere from zero to 
greater than 3 x 10 cm" . Also, the temperature range is T cff = 1.2±g;| x 10 6 K (90% confidence) where 
T e jf is observed at infinity. The large error ranges for T g and Nh prompt us to try to reduce parameter 
uncertainties by better constraining Nh, within the context of this same atmospheric emission model, by 
considering an independent analysis of the same direction. 



The value of Nh obtained from the LAB Survey of Galactic HI (jKalberla et al.l i2005) for this direction 
in the Galaxy, 1.32 xlO 21 cm" 2 , is well within the wide range of acceptable column densities obtained in 
the above model fit. If we now fix Nh at this value in the same absorbed neutron star atmospheric model 
we obtain a new fit with reduced x 2 of 1.000 with 18 de grees of freedom and more precise error ranges: 
T c ff = 8.lj^ 4 x 10 5 K (90% confidence) where T e g is again observed at infinity. The model normalization 
is 1.8lif;?g x 10~ 2 (90% confidence). The derived unabsorbed X-ray flux in the 0.5 to 3 keV range is 
3.1^0 4 x 10~ 14 erg cm" 2 s" 1 (90% confidence). The model normalization gives an indication of the fraction 
of the neutron star surface that is emitting and amounts to a total of ~ 23 km 2 in our simple model, less 
than the entire neutron star surface area. We note that after accounting for the observed background, we 
are working with approximately 300 observed source counts and given this small number and the restricted 
energy range we cannot set strong limits on the column density to the source nor can we investigate the 
possibility of a non-thermal component above 2 keV in the X-ray spectrum. Thus, while it is very likely 
that this X-ray source is in fact the pulsar PSR J2302+4442, longer duration X-ray observations with XMM- 
Newton or Chandra are required to more precisely determine its atmospheric parameters and search for 
possible X-ray pulsations. 



351 4.3. Gamma-ray analysis 

352 The gamma-ray analysis of PSR J2302+4442 was similar to that of PSR J2017+0603 (see Section GDjJ ■ 

353 Figure [2] shows light curves of PSR J2302+4442 in radio and gamma rays. For events within 0.8° of the 

354 MSP the -ff-test parameter is 415.8, also corresponding to a pulsation significance well above 10cr. The 



http: / /heasarc. gsfc.nasa.gov/docs/xanadu/xspcc/ 
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01 23456789 10 

Fig. 4.— The combined EPIC-MOS1 and -MOS2 image of the field of the pulsar PSR J2302+4442, based 
on 24.9 and 25.1 ks exposures, respectively, and smoothed by 3 pixel widths (~ 3.3"). The color scale 
represents counts per pixel. The position of the pulsar is shown by the green cross and is indistinguishable 
to the accuracy of the X-ray image from the position of the X-ray source we call XMMUJ230247+444219. 
This source, and the source labeled A, were both detected by the Swift XRT in its exploration of this field 
(see text) but the two other labeled sources (B & C) apparently were not detected by the XRT. 
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maximum of the radio profile at 1.4 GHz is at phase $ r = 0.960, under the same convention for the absolute 
phasing as described in Section 15751 We checked whether the structure between phase 0.25 and 0.4 comprises 
one or two gamma-ray peaks by plotting light curves with 10, 20, 30, 50 and 100 counts in each bin. We 
found that a sharp peak at phase ~ 0.31 is clearly observed, whereas the possible component at phase ~ 
0.35 is not significant with the present dataset. We fitted the sharp structure at phase ~ 0.31 as well as 
the second gamma-ray peak with Lorentzian functions above constant background. The peak positions and 
FWHM, as well as the radio-to-gamma-ray lag and gamma-ray peak sepa ration are listed in Table El As for 
PSR J2017+0603, the 6 and A values follow the trend already noted bv lAbdo et al.1 (|2010fl ) for previously 
detected gamma-ray pulsars with known radio emission. However, we note in the case of PSR J2302+4442 
an alignment between the radio interpulse at phase ~ 0.65 in Figure [2j and the second gamma-ray peak, 
indicating interesting frequency-dependence of the emission regions, if the radio and the gamma-ray emission 
features are indeed of common origin in the magnetosphere. 

The gamma-ray spectral parameters for PSR J2302+4442 obtained from a fit with j3 = 1 are listed in 
Table[2] and Figure [5] shows the corresponding energy spectrum. In this case, spectral parameters of sources 
within 6° from the pulsar were left free in the fit. The simple power-law model without cutoff is rejected at 
the 9cr level. A spectral fit with the (5 parameter in Equation (TTJ) left free gave f3 = 2.4 ± 0.7. This value 
formally departs from the = 1 assumption; however, we found that there is no statistical improvement of 
the fit compared to the simple exponentially cutoff power-law fit with the current data. As can be seen in 
Figure [5l the best- fit models with j3 = 1 and (3 left free agree well except at the lowest and highest energies, 
where only upper limits could be measured. More data are thus needed to discriminate between the two 
mo dels. Spectral parameters m easured for (3 — 1 are again similar to those of gamma-ray MSPs observed so 
far (jAbdo et al.ll2009ai l2010bh . Finally, the energy flux liste d in Table [2] i s consi stent with that of the 1FGL 



Catalog source J2302.8+4443 measured above 0.1 GeV by lAbdo et all (|2010dl of (4.8 ± 0.4) xlO -11 erg 



We therefore conclude that 1FGL J2302. 8+4443 is associated with the gamma-ray millisecond 
pulsar PSR J2302+4442. 



Discussion 



5.1. Gamma-ray light curve modeling 
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Several of the MSPs detected by the Fermi LAT in gamma rays have quite complex radio pulses, and 
PSRs J2017+0603 and J2302+4442 are no exception. In cont rast, their respective gam ma-ray light curves are 
quite standard, exhibiting a familiar double-peak structure ([Abdo et al.ll2009al l2010ft) . The gamma-ray and 
radio pulse shapes and relative lags motivated light curve modeling using standard outer magnetospheric 
pulsar models commonly employed to describe the light curves of younger pulsars and which have been 
successful in modeling earlier detected gamma-ray MSPs (jVenter et al.1 120091 ) . In such models the gamma- 
ray emission originates in gaps along the last open magnetic field lines, with emission from trailing field 
lines accumulating around a particular o bserver phase leading to intense peaks or "caustics" , due to special 
relativistic effects (iDvks fc Rudakll2003h. In the outer gap (OG) model, two caustics originate from one 
magnetic pole (e.g.. lRomani fc Yadigaro glu 1995), while caustics from both magnetic poles are visible in the 
case of the two-pole caustic (TPC) mode l One may additiona lly consider a pair-starved polar cap (PSPC) 
model ( Muslimov fe Harding 2004 . 20091 : Harding et al. 2005 ) where the combination of perpendicular B- 
field strength and gamma-ray energies of the radiated photons are too low to lead to significant amounts of 
electron-positron pairs close to the stellar surface. In this case, the magnetosphere is "pair-starved" and no 
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Fig. 5. — Same as Figure [31 for PSR J2302+4442. See text for details on the spectral analysis. 
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pair formation front is established, so that the primaries continue to accelerate along the B-field lines and 
emit curvature gamma-ray radiation up to near the light cylinder. The non-zero lags bet ween the gamma-ra y 
and radio pulses led us to model the radio using a phenomcnological model proposed bv lStorv et al. 
where one assumes that the radio emission originates in a cone beam centered on the magnetic dipole axis 
at a single altitude. Different combinations of inclination and observer angles a and £ will result in zero, one 
or two radio peaks from each pole, along with different gamma-ray profile shapes, depending on how close 
an observer's line-of-sight sweeps with respect to the magnetic axis. 

We have used a Markov chain Monte Carlo (MCMC) maximum likelihood fitting technique to jointly 
model the gamma-ray and radio p ulse profiles in order to statistically pick the best-fit emission model 
geometry (details will be described in Johnson et al.ll2011). Additio nally, we have generated simulations with 
1° resolution in a as opposed to the 5° used in I Venter et al. ( 2009 ) and included the Lorentz transformation 
of the magnetic field from the inertial observer's fra me to the co-rotating frame which was missing in previous 
studies and advocated by iBai fc Spitkovskvl (|2010l ) as necessary for self-consistency. An MCMC technique 
involves taking random steps in parameter space, evaluating the likelihood at that step, and accepting the 
step based on th e likelihood ratio with the previous step. In particular, we use a Metropolis-Hastings method 
(|Hastingslll9"70h to update the parameter state, accepting steps if the likelihood at the new step is greater 
than the previous step or if the ratio is greater than a random number £ [0, 1). For eac h model fit we verify 
that our MCMC has converged using the method proposed bv lGelman fe Rubinl (|1992f ) . 



The gamma-ray light curves are fit using Poisson likelihood and the radio profiles using a x 2 statistic. 
In order to balance the contributions from the radio and gamma-ray data, and in particular to balance the 
high statistical precision of the radio data against our simple cone-beam model, we have used a relative error 
for the radio data equal to the average gamma-ray relative uncertainty in the on-peak region times the radio 
maximum. It is important to note that the choice of uncertainty for the radio profile can strongly affect the 
best-fit results. A smaller uncertainty will decrease the overall likelihood, which can in some cases lead to a 
different best-fit geometry favoring the radio light curve. For both MSPs we have taken the gamma-ray on- 
peak interval to be <f> 6 [0.25,0.75]. Our geometric models assume constant-emissivity gamma-ray emission 
extending from the stellar surface in the TPC model, while the minimum radius is set to the radius of the null 
charge surface (which depends on magnetic azimuth and co-latitude) in the OG model. For all simulations 
we have used a maximum emission altitude for the gamma rays of 1.2 RlCi where Rlc = cP/(2tt), with 
the added caveat that the emission not go beyond a cylindrical radius equal to 0.95 Rlc- We found that 
the likelihood surfaces are very multi-modal which can lead to a low acceptance rate and an incomplete 
explo ration of the parameter space; therefore, we have implemented simulated tempering (IMarinari fe Parisi 
1992) with small-world chain steps ( Guan et al. 20061) in a and C- The MCMC parameter space includes a, 



C (both with 1° resolution), gap width w (with a resolution of 0.05, normalized to the polar cap radius), and 
phase-shift, which accounts for the fact that the definitions of phase zero are different between the data and 
our models. Our MCMC is implemented in python using the scipy moduk@ and the li ght curve fitting for 
each step is done using the scipy. optimize. fmin_Lfbgs-b multi-variate, bound optimizer (|Zhu et al.lll997l) . 



In order to match the data with our simulations we re-binned both the gamma-ray and radio data to 
60 bins, see Figures [6] and [71 This has the effect of smoothing out very fine scale variations in the radio 
profile, but as we discuss below our radio profile simulations arc not refined enough to reproduce these 
structures and thus fitting to the 60 bin radio profiles is sufficient to reproduce the general features, namely 
the gamma-to-radio lag. For PSR J2017+0603 we find best-fit solutions of a = 16° and £ = 68° with an 



7 See |http : / /d ocs . scipy. org/ doc/] for documentation 
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449 



infinitely thin gap for a TPC model and a — 17° and £ = 68° with an infinitely thin gap for an OG model. 
For PSR J2302+4442 we find best-fit solutions of a = 58° and C = 46° with infinitely thin gap for a TPC 
model and a — 63° and £ = 39° with infinitely thin gap for an OG model. When we find best-fit models 
with infinitely thin gap widths for both pulsars we do not think this represents the truth as a zero-width 
gap is unphysical; rather, we take this to mean that the best gap width is somewhere between and 0.05 
and the best-fit value of is chosen only as a result of the resolution of our simulations. Note also that we 
have not yet calibrated the fitting procedure to address the significance of differences in — log(likelihood) so 
we cannot be more quantitative in discussing the preference of one model over another. However, for both 
MSPs differences in — log(likelihood) were close to 0, meaning that neither of TPC and OG geometries are 
preferred. 

Neither of the model fits for PSR J20 17+0603 are able to produce a wide enough first gamma-ray peak 
but both produce the correct peak separation. Also, the model fits cannot reproduce all the features observed 

450 in the radio profile. However, the best-fit geometries are able to produce radio-to-gamma-ray lags close to 

451 what is observed. The situation is similar for PSR J2302+4442, with both models matching the sharp second 

452 gamma-ray peak but neither is able to produce a strong enough first peak. The TPC model implies two 

453 small peaks near phase 0.3 for slightly different values of a and close the to best-fit values. Tests have 

454 shown that lowering the maximum emission altitude can affect the prominence of these two peaks, which 

455 suggests that more investigation is merited in this parameter. With more data the significance, or not, of 

456 this two-peaked structure will serve as a further discriminator between the models. Neither best- fit geometry 

457 produces produces two radio peaks with the correct spacing. The TPC geometry does predict two closely 
45a spaced radio peaks while the OG geometry approximately matches the radio peak near 0.15 in phase. 

459 For both MSPs, it is of interest to note that geometries with a and £ both near 20° produce two 

460 radio peaks with approximately correct spacing but the resultant gamma-ray TPC light curves are similar 

461 to square waves while the gamma-ray emission in OG models is missed entirely. Clearly, our simple radio 

462 model does not adequately reproduce the data. Both MSPs have at least three components in their radio 

463 profiles, while the model can only produce zero, one, or two peaks from each magnetic pole. This points to 

464 more complex radio emission geometries, with radio emission from both magnetic poles visib l e, and likely 

465 that emission may occur higher up in the magnetosphere as has been suggested bv lRavi et al. ( 2010l ). 



We also fit both MSPs with the PSPC model, though this is not as successful at producing sharp 
gamma-ray peaks. For both MSPs the fits predict a ~ 70° and ( ~ 80° which suggest that we would see 
radio emission from both magnetic poles. The gamma-ray PSPC models are able to reproduce the second, 
sharp peak for each MSP but have trouble matching the first peak properly. The best-fit geometries result 
in more complex radio profiles but are still not able to match all of the observed features. For both MSPs 
the PSPC models are disfavored by the likelihood when compared to the TPC and OG fits. Our modeling 
and fit results also show that there is still much to be learned about the radio beam structure. 



473 5.2. Gamma- ray efficiencies 

474 One can derive the total gamma-ray luminosity above 0.1 GeV and the efficiency of conversion of spin- 

475 down energy into gamma rays with the following expressions: 



L-y = AirfnGd 2 , (2) 
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Fig. 6. — Top: gamma-ray data and modeled light curves for PSR J2017+0603 with 60 bins per rotation. 
Bottom: Nangay 1.4 GHz radio profile and modeled light curves. Modeled light curves were made using 
a = 16°, £ = 68° and an infinitely thin gap for the TPC model, and a — 17°, C = 68° and an infinitely thin 
gap for the OG geometry. See Section UTTl for emission altitude extents. 



-20- 




Fig. 7 — Same as Figure El for PSR J2302+4442. Modeled light curves were made using a = 58°, C = 46° 
and an infinitely thin gap for the TPC emission geometry, and a = 63°, £ = 39° and an infinitely thin gap 
for the OG model. 
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T) = L./E. (3) 

In these expressions, d and E are the pulsar distance and spin-down energy, /q is the correction factor 
depending on the viewing geometry defined above, and G is the energy flux measured above 0.1 GeV. Table 
[5] lists Lry and 7/ values under the assumption that fa = 1, and using the pulsar distances inferred from 
the NE2001 model (see Table [T|). For both pulsars, gamma-ray efficiencies are found to be suspiciously 
large, and even greater than 100% in the case of PSR J2302+4442, which is unphysical. Overestimated 
/n factors and distances are plausible explanations for the large efficiency values. The best-fit TPC and 
OG emission geometries discussed in Section 15.11 predict geometrical correction factors of 0.48 and 0.30 
for PSR J2017+0603, leading to realistic gamma-ray efficiencies of 0.39 and 0.24, respectively. However, 
fn factors calculated under TPC and OG geometries for PSR J2302+4442 are 0.95 and 0.97, leading to 
gamma-ray efficiencies greater than 1.6. The distance inferred from the NE2001 model is therefore likely 
overestimated, or the model is incorrect. In addition, proper motions could make the apparent spin-down 
energy loss rates E larger than the intrinsic values because of the Shklovskii effe ct, thereby increasing gamm a- 



ray efficiencies. The average efficiency of gamma-ray MSPs observed so far (jAbdo et al.ll2009al |2010bJ) is 
~ 10% (we excluded PSR J1614— 2230, which also has an unphysical gamma-ray efficiency of 100% with 
the NE2001 distance). Assuming an efficiency of 10% for PSR J2302+4442, we find that the distance has 
to be smaller by a factor of 4, which would place the pulsar at d < 300 pc. Note however that the X- 
ray energy flux Gx of ~ 3.1 x 10~ 14 erg cm -2 s _1 measured between 0.5 and 3 keV leads to an X-ray 
efficiency of AitGxd 2 / E ~ 1.4 x 10~ 3 if we assume the NE2001 distance of 1.18 kpc, while it decreases to 
~ 9 x 10~ 5 with a distance of 3 pc. The former efficiency is very close to the 10 -3 value empirically 



predicted by iBecker fc Truemperl (jl997l ) at these energies. The X-ray analysis therefore does not support 



such an important reduction of the distance. If the pulsar distance is indeed that small, a timing parallax 
7r = d ( k p C ) ^ 3.3 mas should be measurable with accumulated radio timing observations. This parallax could 
also be measured via the VLBI measurements being undertaken for all Fermi pulsar^f). 



6. Conclusions 

In a search for radio pulsations at the position of Fermi 1FGL catalog sources with the Nancay radio 
telescope, we discovered two millisecond pulsars, PSRs J2017+0603 and J2302+4442, both orbiting low- 
mass companion stars. Both pulsars were found to emit pulsed gamma-ray emission, indicating that they 
are associated with the previously unidentified gamma-ray sources. The gamma-ray light curves and spectral 
properties of the two MSPs are reminiscent of those of other gamma-ray MSPs observed previously. 

Prior to Fermi, error boxes of unidentified gamma-ray sources were much larger than radio telescope 
beams, making searches for pulsars dif ficult, as multiple poin tings were required to cover the gamma-ray 



source contour entirely (see for example IChampion et al.ll2005l ). Unassociated Fermi LAT sources are typi- 



cally localized to within 10 arcminutes, which is comparable to radio beam sizes and therefore makes radio 
pulsation searches easier and more efficient. With its improved localization accuracy and its homogeneous 
coverage of the gamma-ray sky, the Fermi LAT is therefore revealing the population of energetic pulsars 
and millisecond pulsars, providing a complementary view of the Galactic population of pulsars, which has 
mostly been studied at radio wavelengths up to now. 



Cycle 3 Fermi Guest Investigator proposal: S. Chatterjee et al. 
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Table 1: Parameters for PSRs J2017+0603 and J2302+4442. See Sections EQ and E0 for details on the 

measurement of these parameters. Numbers in parentheses are the nominal la tempo2 uncertainties in the 
least-significant digits quoted. 

Parameter PSR J2017+0603 

Right ascension (J2000) 20:17:22.7044(1) 

Declination (J2000) 06:03:05.569(4) 

Rotational period, P (ms) 2.896215815562(2) 

Period derivative, P (10~ 21 ) 8.3(1) 

Epoch of ephemeris, T (MJD) 55000 

Dispersion measure, DM (cm -3 pc) 23.918(3) 

Orbital period, P b (d) 2.198481129(6) 

Projected semi-major axis, x (It s) 2.1929239(7) 

Epoch of ascending node, T asc (MJD) 55202.5321589(3) 

esinu; 0.0000023(6) 

ecosu -0.00000046(6) 

Span of timing data (MJD) 54714 — 55342 

Number of TO As 71 

RMS of TOA residuals (jits) 3.23 

Units TDB 

Solar system ephemeris model DE405 

Flux density at 1.4 GHz, g 14 oo (mjy) 0.5(2) 

Derived parameters 

Orbital eccentricity, e 0.000005(2) 

Mass function, / (M ) 0.002342653(2) 

Minimum companion mass, m c (Mq) > 0.18 

Galactic longitude, I (°) 48.62 

Galactic latitude, 6 (°) —16.03 

Distance inferred from the NE2001 model, d (kpc) 1.56 ± 0.16 

Spin-down luminosity, E (10 33 erg s _1 ) 13.43 

Characteristic age, r (10 9 yr) 5.55 

Surface magnetic field strength, B s (10 s G) 1.57 

Magnetic field strength at the light cylinder, Blc (10 4 G) 5.86 



PSR J2302+4442 
23:02:46.9796(7) 
+44:42:22.090(5) 
5.192324646411(7) 
13.3(5) 
55000 
13.762(6) 
125.935292(3) 
51.429942(3) 
55096.517187(3) 
-0.00023537(6) 
-0.00044485(6) 
54712 — 55342 
130 
6.46 
TDB 
DE405 
1-2(4) 

0.0005033(2) 
0.009209497(1) 
> 0.30 
103.40 
-14.00 

i lo+O.lt) 

J - 1H -o.a3 
3.74 

6.20 

2.66 

1.73 
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Table 2: Light curve and spectral parameters of PSRs J20 17+0603 and J2302+4442 in gamma rays, fixing 
j3 = 1 in Equation (JlJ. See Sections 13.31 and 14.31 for details on the measurement of these parameters. Peak 
positions, widths and separations are given in phase units, between and 1. 

Parameter PSR J2017+0603 PSR J2302+4442 

First peak position, $i 0.348 ± 0.009 0.310 ± 0.021 

First peak full width at half-maximum, FWHMi 0.248 ± 0.054 0.033 ± 0.013 

Second peak position, $ 2 0.636 ± 0.005 0.629 ± 0.003 

Second peak full width at half-maximum, FWHM 2 0.050 ± 0.013 0.037 ± 0.006 

Radio-to-gamma-ray lag, 5 0.225 ± 0.009 ± 0.002 0.350 ± 0.021 ± 0.002 

Gamma-ray peak separation, A 0.288 ± 0.010 0.320 ± 0.021 



Spectral index, F 


1.00 


± 


0.16 


± 


0.16 


1.25 ± 0.13 


± 0.14 


Cutoff energy, E c (GeV) 


3.12 


± 


0.57 


± 


0.75 


2.97 ± 0.51 


± 0.54 


Photon flux, F (> 0.1 GeV) (10~ 8 cm' 2 s' 1 ) 


2.21 


± 


0.31 


± 


0.11 


3.34 ± 0.38 


± 0.20 


Energy flux, G (> 0.1 GeV) (10' 11 erg cm" 2 s" 1 ) 


3.71 


± 


0.24 


± 


0.19 


3.94 ± 0.22 


± 0.10 


Luminosity, L 7 / fn (10 33 erg s _1 ) 


10.79 


± 


1.72 


± 


1.66 


fi c: 7 +0.87 
D - ' -1.85 


+0.80 

-1.82 


Efficiency, rj / fn 


0.80 


± 


0.13 


± 


0.12 


i 7 r +0.23 
1 - ' -0.49 


+0.21 
-0.48 



